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Gamow- Teller strength distributions, /?-decay half-lives, and /3-delayed neutron emission are inves- 
tigated in neutron-rich Zr and Mo isotopes within a deformed quasiparticle random-phase approx- 
imation. The approach is based on a self-consistent Skyrme Hartree-Fock mean field with pairing 
correlations and residual separable particle-hole and particle-particle forces. Comparison with re- 
cent measurements of half-lives stresses the important role that nuclear deformation plays in the 
description of /3-decay properties in this mass region. 
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I. INTRODUCTION 

The rapid neutron-capture process (r process) is con- 
sidered as the main nucleosynthesis mechanism respon- 
sible for the production of heavy neutron-rich nuclei and 
for the existence of about half of the nuclei heavier than 
iron 0, Although the astrophysical sites for this 
process are still controversial, it takes place in scenar- 
ios characterized by very high neutron densities. The 
path that nucleosynthesis follows involves neutron-rich 
isotopes, which can be far away from the valley of in- 
stability. The most relevant nuclear properties to de- 
scribe the r process are the nuclear masses and the /3- 
decay properties [H, H[ — namely, the /3-decay half-lives 
(T1/2) and the /3-delayed neutron-emission probabilities 
(P„). Nuclear masses define the possible r-process paths 
near the neutron drip-lines. The T\ii values of r-process 
waiting-point nuclei determine the pre freeze-out isotopic 
abundances and the speed of the process towards heav- 
ier elements, as well as the r-process time scale. The P n 
values of r-process isobaric nuclei define the decay path 
towards stability following the freeze-out, and provide a 
source of late-time neutrons. 

A reliable nuclear physics description of the properties 
of the extremely neutron-rich nuclei along the r-process 
path is needed to interpret the astrophysical observations 
and to model and simulate properly the r process. The 
quality of the nucleosynthesis modeling is directly af- 
fected by the quality of the nuclear structure input. Un- 
fortunately, most of the nuclear properties of relevance 
for the r process are experimentally unknown, although 
much effort is being done recently, and therefore theo- 
retical predictions must be considered. Such calculations 
are particularly challenging in the very exotic regions of 
interest, since they involve extrapolations using well es- 
tablished nuclear-structure models that have been prop- 
erly tuned to account mostly for the properties of nuclei 
in the valley of stability. In particular, the shell structure 
of neutron-rich drip-line nuclei is still unknown to a large 
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extent. Significant isospin dependence of shell effects in 
medium- mass and heavy nuclei has been predicted [343 ■ 
It has been found that the shell gaps dramatically de- 
crease near the neutron drip-lines because of continuum 
effects and a quenching of shell effects is apparent. 

As a matter of fact, nuclear structure properties of 
nuclei far from stability where no experiments exist for 
direct comparison, can be tested by exploring their in- 
fluence on the solar r-process abundance patterns. As 
an example, the agreement with the observed r-process 
abundances in the A ~ 120 mass region is manifestly im- 
proved 0-G3 

when using nuclear structure models that 
include a shell quenching effect at N = 82 0, [ll[ . 

In this work we focus our attention on the mass re- 
gion of neutron-rich A ~ 100 — 110 nuclei, which is of 
great interest for the astrophysical r process. In addition, 
neutron-rich isotopes in this mass region are known [l2| 
to be interesting examples where the equilibrium shape 
of the nucleus is rapidly changing and shape coexistence 
is present with competing prolate, oblate, and spherical 
shapes at close energies (see e.g. 
view) . 

In a recent publication )14j, the /3-decay properties of 
some neutron-rich Zr and Mo isotopes were measured for 
the first time. The data were interpreted in terms of 
the quasiparticle random-phase approximation (QRPA) 
[l5l4l8| ] , using nuclear shapes and nuclear masses derived 
from the finite-range droplet model (FRDM) Q1] and 
the latest version of the finite-range liquid-drop model 
(FRLDM) 20], which also includes triaxial deforma- 
tion. QRPA calculations for neutron-rich nuclei have also 
been performed within different approaches, such as the 
Hartree-Fock-Bogoliubov (HFB) [Hj], continuum QRPA 
with either the extended Thomas-Fermi plus Strutinsky 
integral (ETFSI) method [22j or based on density func- 
tional f2 3| . [2 4l. and the relativistic mean field (RMF) 
approach |25|. iust to mention some of the recent publica- 
tions, all of them for spherical nuclei. However, the mass 
region of concern here requires nuclear deformation as 
a relevant degree of freedom to characterize the nuclear 
structure involved in the calculation of the /3-strength 
functions. The deformed QRPA formalism has been de- 
veloped in Refs. [TB-fU HSH3 1 where phenomenological 
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mean fields based on Nilsson or Woods-Saxon potentials 
were used as a starting basis. In this work we investigate 
the decay properties of neutron-rich even-even Zr and 
Mo isotopes within a deformed self-consistent Hartree- 
Fock (HF) mean field formalism with Skyrme interactions 
and pairing correlations in BCS approximation. Residual 
spin-isospin interactions are also included in the particle- 
hole a nd p article-particle channels and are treated in 

qrpa H2HJ. 

The paper is organized as follows. In Sec. [TT| a brief 
review of the theoretical formalism is presented. Sec. 
IIIII contains the results obtained within this approach for 
the potential energy curves, Gamow- Teller (GT) strength 
distributions, and /3-decay half-lives. Sec. IV summarizes 
the main conclusions. 



II. THEORETICAL FORMALISM 

In this section we show briefly the theoretical frame- 
work used in this paper to describe the /3-decay prop- 
erties in Zr and Mo neutron-rich isotopes. More details 
of the formalism can be found in Refs. [29, 30]. The 
method consists of a self-consistent formalism based on a 
deformed Hartree-Fock mean field obtained with Skyrme 
interactions, including pairing correlations. The single- 
particle energies, wave functions, and occupation proba- 
bilities are generated from this mean field. In this work 
we have chosen the Skyrme force SLy4 [3l| as a repre- 
sentative of the Skyrme forces. This particular force in- 
cludes some selected properties of unstable nuclei in the 
adjusting procedure of the parameters. It is one of the 
more successful Skyrme forces and has been extensively 
studied in the last years. 

The solution of the HF equation is found by using the 
formalism developed in Ref. 32] , assuming time reversal 
and axial symmetry. The single-particle wave functions 
are expanded in terms of the eigenstates of an axially 
symmetric harmonic oscillator in cylindrical coordinates, 
using twelve major shells. The method also includes pair- 
ing between like nucleons in BCS approximation with 
fixed gap parameters for protons and neutrons, which are 
determined phenomenologically from the odd-even mass 
differences through a symmetric five term formula involv- 
ing the experimental binding energies (33| when avail- 
able. In those cases where experimental information for 
masses is still not available, we have used the same pair- 
ing gaps as the closer isotopes measured. The pairing 
gaps for protons (A p ) and neutrons (A„) obtained in 
this way are roughly around 1 MeV. The correspond- 
ing pairing strengths G p and G n calculated from the gap 
equation depend sensitively on the mass region, single- 
particle spectrum, and active window for pairing. For 
typical values of the cutoffs of about 5 MeV around the 
Fermi level, one obtains G p ~ 0.25 MeV and G n ~ 0.30 
MeV. It is worth noticing that, although the BCS formal- 
ism leads to an unphysical neutron gas surrounding the 
nucleus near the drip line, the approximation is still valid 



in the region considered here, where the pairing gaps are 
still much lower than the Fermi energies. 

The potential energy curves (PEC) are analyzed as a 
function of the quadrupole deformation /3, 

written in terms of the mass quadrupole moment Qo and 
the mean square radius (r 2 ). For that purpose, con- 
strained HF calculations are performed with a quadratic 
constraint [34]. The HF energy is minimized under the 
constraint of keeping fixed the nuclear deformation. Cal- 
culations for GT strengths are performed subsequently 
for the equilibrium shapes of each nucleus, that is, for 
the solutions, in general deformed, for which minima are 
obtained in the energy curves. Since decays connect- 
ing different shapes are disfavored, similar shapes are as- 
sumed for the ground state of the parent nucleus and for 
all populated states in the daughter nucleus. The valid- 
ity of this assumption was discussed for example in Refs. 

pint In 

our particular case, for SLy4 and neutron- 
rich Zr and Mo isotopes, the ground-state deformation of 
the even-even parents (Zr,Mo) and of the corresponding 
/3-decay odd-odd daughters (Nb,Tc) are practically the 
same, as it can be seen in Ref. [35j. 

To describe GT transitions, a spin-isospin residual in- 
teraction is added to the mean field and treated in a de- 
formed proton- neutron QRPA [H-[l3 • This inter- 
action contains two parts, particle-hole (ph) and particle- 
particle (pp). The interaction in the ph channel is re- 
sponsible for the position and structure of the GT res- 
onance [26|, [5^, [301 and it can be derived consistently 
from the same Skyrme interaction used to generate the 
mean field, through the second derivatives of the energy 
density functional with respect to the one-body densi- 
ties. The ph residual interaction is finally expressed in a 
separable form by averaging the Landau-Migdal result- 
ing force over the nuclear volume, as explained in Refs. 
[2^, ■ The coupling strength is given by Xgt = 0.15 
MeV. The pp part is a neutron-proton pairing force in the 
J v = 1 + coupling channel, which is also introduced as a 
separable force [U [H[ . The strength of the pp residual 
interaction in this theoretical approach is not derived sclf- 
consistently from the SLy4 force used to obtain the mean 
field basis, but nevertheless it has been fixed in accor- 
dance to it. This strength is usually fitted to reproduce 
globally the experimental half-lives. Various attempts 
have been done in the past to fix this strength, arriv- 
ing to expressions such as k P q T = 0.58/j4 ' 7 in Ref. [26| . 
which depend on the model used to describe the mean 
field, Nilsson model in the above reference. Our work in 
the past (see Refs. [H, [13] and references therein) , based 
on the Skyrme force SLy4 leads us to consider the value 
k P q T = 0.03 MeV as a reasonable choice in this mass 
region, improving the agreement with the experimental 
half-lives. 
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The proton-neutron QRPA phonon operator for GT 
excitations in even-even nuclei is written as 



r+ 

UK 



(2) 



where a + (a) are quasiparticle creation (annihilation) op- 
erators, lok are the QRPA excitation energies with re- 
spect to the ground state of the parent nucleus, and 
X*£*f ,Y"k the forward and backward amplitudes, re- 
spectively. For even-even nuclei the allowed GT tran- 
sition amplitudes in the intrinsic frame connecting the 
QRPA ground state |0) (T utc |0) =0) to one-phonon 
states \u> K ) (T+ k |0) = \uik)), are given by 



where 



with 



Sun 
K i 



(3) 

(4) 
(5) 

(6) 



v's are occupation amplitudes (u 2 = 1 — v 2 ) and spin 
matrix elements connecting neutron and proton states 
with spin operators 



SUIT 
K 



(v\a K \n). (7) 

The GT strength B (Aj (GT ± ) in the laboratory system 
for a transition 7ji£j(0 + 0) — > I fK f(l + K) can be ob- 
tained in terms of the intrinsic amplitudes in Eq. ([3]) 

as 

B^GT*) = J2 [(^K=o h^O) 2 *(w K =o - w) 

U) K 

+2{u) K=l \ai^\Q) 2 6{u K= i-u^ ,(8) 

in [g A /4ir] units. To obtain this expression, the initial 
and final states in the laboratory frame have been ex- 
pressed in terms of the intrinsic states using the Bohr- 
Mottelson factorization (38|. 

The excitation energy E ex referred to the ground state 
of the odd-odd daughter nucleus is obtained by sub- 
tracting the lowest two-quasiparticle energy E from the 
calculated uj energy in the QRPA calculation, E ex = 
ojqrpa — Eq, where Eq = (E n + E p )o is the sum of the 
lowest quasiparticle energies for neutrons and protons. 
The GT strength B(GT) will be plotted later versus E ex 
in Figs. 1H El and El 

The /3-decay half-life is obtained by summing all the 
allowed transition strengths to states in the daughter nu- 
cleus with excitation energies lying below the correspond- 
ing Q-energy, and weighted with the phase space factors 
f(Z,Q fj -E ex ), 



1 l/2 ~ 



(9A/gv) eS 



D 



E 



f(Z,Q p -E ex )B(GT,E e 



(9) 

with D = 6200 s and (gA/gv)cS = 0.77(gA/gv)hcc, 
where 0.77 is a standard quenching factor that takes into 
account in an effective way all the correlations which 
are not properly considered in the present approach. The 
bare results can be recovered by scaling the results in this 
paper for B(GT) and T±/2 with the square of this quench- 
ing factor. The Qp- energy is given by 



Qp- = M(A,Z)-M(A,Z+l)-m e 
= BE(A, Z) - BE(A, Z + 1) + 1 



(10) 



written in terms of the nuclear masses M(A, Z) or nuclear 
binding energies BE(A, Z) and the neutron (m n ), proton 
(jn p ), and electron (m e ) masses. 

The Fermi integral f(Z, Qp — E ex ) is computed numer- 
ically for each value of the energy inclu ding screening and 
finite size effects, as explained in Ref. pOj . 

f p (Z,W Q )= J pW(W Q - W) 2 X ± (Z,W)AW , (11) 
with 



\±(Z, W) = 2(1 + 7 )(2pfl)-^-7) e Tn> ffl±^ , 

7 (12) 

where 7 = yl — (aZ) 2 ; y — aZW/p ; a is the fine struc- 
ture constant and R the nuclear radius. W is the total 
energy of the ft particle, Wp is the total energy available 
in m e c 2 units, and p = \JW 2 — 1 is the momentum in 
m e c units. 

This function weights differently the strength B(GT) 
depending on the excitation energy. As a general rule 
f(Z, Qj3~E ex ) increases with the energy of the /3-particle 
and therefore the strength located at low excitation en- 
ergies contribute more importantly to the half-life. 

The probability for /3-delayed neutron emission is given 

by 



J2 f(Z,Q l3 -E ex )B(GT,E ex ) 



Pn 



S n <E cx «. 



J2 f(Z,Q l3 ~E ex )B(GT,E ex ) 

0<E^<Qf, 



(13) 



where the sums extend to all the excitation energies in 
the daughter nuclei in the indicated ranges. S n is the 
one-neutron separation energy in the daughter nucleus. 
In this expression it is assumed that all the decays to 
energies above S n in the daughter nuclei always lead to 
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delayed neutron emission and then, 7-decay from neu- 
tron unbound levels is neglected. According to Eq. (fl"3"|). 
P n is mostly sensitive to the strength located at energies 
around S n , thus providing a structure probe complemen- 
tary to Ti/2- 



III. RESULTS 

In this section we start by showing the results obtained 
for the potential energy curves in the isotopes under 
study. Then, we calculate the energy distribution of the 
GT strength corresponding to the local minima in the 
potential energy curves. After showing the predictions 
of various mass models to the Qp and S n values for the 
more unstable isotopes, where no data on these quanti- 
ties are available, we calculate the /3-decay half-lives and 
discuss their dependence on the deformation. 

In previous works we have studied the 

sensitivity of the GT strength distributions to the various 
ingredients contributing to the deformed QRPA-like cal- 
culations, namely to the nucleon-nucleon effective force, 
to pairing correlations, and to residual interactions. We 
found different sensitivities to them. In this work, all of 
these ingredients have been fixed to the most reasonable 
choices found previously and mentioned above, including 
the quenching factor. Here, we mainly discuss effects of 
deformation, keeping in mind that the method provides 
the self-consistent deformations as well. 
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A. Potential Energy Curves 

In Fig. [2[2]) we can see the potential energy curves 
for the even-even i°o-u°Zr ( 104 - 114 Mo) isotopes. We 
show the energies relative to that of the ground state 
plotted as a function of the quadrupole deformation /3. 
They are obtained from constrained HF+BCS calcula- 
tions with the Skyrme forces SG2 [45] and SLy4 [Hj]. We 
observe that both forces produce very similar results. In 
Fig. [T]we see that the Zr isotopes exhibit in all cases two 
well developed minima. The ground states are located 
in the prolate sector at positive values of (3 ps 0.4. We 
can also see oblate minima at higher energies located at 
ft ps —0.2. The two minima are separated by potential 
energy barriers varying from E = 3 MeV, in the light- 
est 100 Zr isotope, up to barriers of the order of 5 MeV 
in heavier isotopes. In the isotopes 108 - 110 Zr a spherical 
local minimum is also developed. 

Similar trends are observed in Fig. [5] for the Mo iso- 
topes. We observe well developed oblate and prolate min- 
ima, which are separated by barriers ranging from 3 MeV 
up to 5 MeV. We get a prolate ground state with an 
oblate minimum very close in energy in the lightest iso- 
tope considered, 104 Mo, a practically degenerate oblate- 
prolate in 106 Mo, and oblate ground states in heavier 
isotopes with quadrupole deformations at (3 ps —0.2 with 
prolate excited states at energies lower than 1 MeV. 



FIG. 1: Potential energy curves for even-even Zr iso- 

topes obtained from constrained HF+BCS calculations with 
the Skyrme forces SG2 and SLy4. 



Again, the heavier isotopes favor the appearance of a 
spherical configuration at very low energies, resulting in 
an emergent triple oblate-spherical-prolate shape coexis- 
tence scenario. 

These results are in qualitative agreement with similar 
ones obtained in this mass region from different theoret- 
ical approaches. As an example of these methods we can 
mention the results obtained in Ref. (4(| , where this mass 
region was studied within a macroscopic-microscopic ap- 
proach based on an energy obtained from a liquid drop or 
a FRLDM modified by a shell correction taken from a de- 
formed Woods-Saxon potential. Zr isotopes from TV = 60 
to N — 72 were predicted to have well-deformed pro- 
late ground states, while Mo isotopes suffered a shape 
transition from prolate shapes in the lighter neutron-rich 
isotopes (N = 62) to oblate shapes in the heavier ones. 
Similarly, the deformations obtained in Ref. [19( from 
the FRDM and a folded- Yukawa single-particle micro- 
scopic model were in the range (3 — 0.36 — 0.38 in the Zr 
isotopes considered in this work and j3 — 0.33 — 0.36 in 
the Mo isotopes, except in the heavier 114 Mo, where an 
oblate shape with j3 ps —0.25 becomes the ground state. 
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FIG. 2: Same as in Fig. [Q but for 104 - 114 Mo isotopes. 



RMF calculations [471 . |48| show ground-state deforma- 
tions in the range of 8 = 0.36 — 0.40 in the Zr isotopes, 
while for Mo isotopes oblate ground states are obtained 
with parameters of deformation between = —0.28 and 
= —0.23, except in the lighter isotope 104 Mo, where a 
prolate ground state — 0.336 is found. Calculations in- 
cluding rotational states in terms of the total Routhian 
surface (TRS), using non-axial Wood-Saxon potentials 
[ill, predicted two coexisting prolate and oblate min- 
ima {(3 k, 0.35 and (3 w -0.2) for ">6-ii6 Zr i soto pes, 
where the prolate ground state becomes oblate beyond 
110 Zr. The same calculations showed oblate ((3 s» —0.22) 
ground-states for N > 68 Mo isotopes. Finally, simi- 
lar results in the sense of competing oblate and prolate 
shapes and emergence of spherical configurations in the 
heavier isotopes are also obtained within the Hartree- 
Fock-Bogoliubov framework with the finite-range effec- 
tive Gogny interaction D1S [5(|. Equilibrium oblate 
(f3 ss —0.2) and prolate (8 ~ 0.4) coexistent deforma- 
tions were found in Ref. [5f| practically at the same en- 
ergy in the Zr isotopes. In Mo isotopes an oblate shape 
(8 « —0.2) is favored energetically with close prolate 
(0 m 0.4) solutions. In both cases, Zr and Mo isotopes, 
a spherical solution lows in energy and becomes almost 



degenerate with the deformed solutions for the heavier 
isotopes 110 Zr and 114 Mo. 

Thus, a consistent theoretical picture emerges, which 
is supported by the still scarce experimental informa- 
tion available. Experimentally, two coexisting deformed 
bands weakly admixed were found in 100 Zr [Bit I52T] from 
an analysis of B(E2) and p(E0) and a two-level mixing 
model analysis. One of these bands is a highly deformed 
prolate yrast band (0 = 0.34), while the other is mod- 
erately deformed (|/3| = 0.16) and weakly mixed to the 
yrast by about 10%. The highly deformed band in 100 Zr 
is nearly identical to the yrast band in 102 Zr. Hill et al. 
[53l ] have also discussed the possibility that the level 
measured for 102 Zr at 895 keV could be the head of a 
band with \0\ w 0.2, similar to the S band of 100 Zr. 

Quadrupole moments were also determined [54| for ro- 
tational bands in 98 ~ 104 Zr isotopes and deformation pa- 
rameters were deduced increasing gradually from = 0.1 
at N = 56 up to 8 = 0.4 at N = 64. More recently [EH, 
large deformations (8 — 0.47(7)) were extracted in 104 Zr 
and in 106 Mo (8 = 0.36(7)) from the half-lives of their 
2+ states. Spectroscopic studies of high-spin states of 
ioo-io 4 2 r and 102_108 Mo have also been performed by 
Hua et al. [56| within the particle-rotor model. Accord- 
ing to these authors, the difference in signature split- 
ting observed for the 5/2 ~ [532] band between the odd 
Zr and Mo isotopes could be attributed to the appear- 
ance of triaxiality in Mo isotopes. As mentioned above, 
the formalism employed in the present study does not in- 
clude non-axial deformation. Such limitation, however, 
has no significant impact in the results discussed here. 
As an example, the inclusion of triaxiality in the last 
version of the FRLDM 20] resulted in a small reduction 
of the 106 ' 108 Mo ground-state energies (of about 250 keV) 
at 7=17.5°, with respect to pure prolate shapes. Simi- 
larly, Xu et al. [49( predict a 7-soft triaxial minimum for 
108 Mo. 



B. Gamow- Teller strength distributions 

In the next figures, we show the results obtained for 
the energy distributions of the GT strength correspond- 
ing to the oblate-prolate-spherical equilibrium shapes for 
which we obtained minima in the potential energy curves 
in Figs. [T]and[2j The results are obtained with the force 
SLy4, using constant pairing gaps extracted from the ex- 
perimental masses (or systematics) and with residual in- 
teractions with the parameters written in Sec. [II] The 
GT strength in (g^/An) units, is plotted versus the ex- 
citation energy of the daughter nucleus and a quenching 
factor 0.77 has been included. 

Figs. [3] and 0] contain the results for Zr and Mo iso- 
topes, respectively. We show the energy distributions 
of the individual GT strengths together with continu- 
ous distributions obtained by folding the strength with 1 
MeV width Breit-Wigner functions. The vertical arrows 
show the Qp and S n energies, taken from experiment (33j 
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FIG. 3: (Color online) QRPA-SLy4 Gamow-Teller strength 
distributions for Zr isotopes as a function of the excitation 
energy in the daughter nucleus. The calculations correspond 
to the various equilibrium configurations found in the PECs 
curves. Qp and S n values are shown by solid and dashed 
vertical arrows, respectively. 



or from the mass formula in Ref. [57| when data are not 
available, as we shall explain later on. 

The main characteristic of these distributions is the ex- 
istence of a GT resonance located at increasing excitation 
energy as the number of neutrons N increases. The to- 
tal GT strength also increases with TV, as it is expected 
to fulfill the Ikeda sum rule. It is worth noticing that 
both oblate and prolate shapes produce quite similar GT 
strength distributions on a global scale. Even the spher- 
ical profiles are quite close to the deformed ones. Never- 
theless, the small differences among the various shapes at 
the low energy tails (below the Qp) of the GT strength 
distributions that can be appreciated because of the log- 
arithmic scale, lead to sizable effects in the /3-decay half- 
lives. In the next figures, Fig. [5]for Zr isotopes and Fig. [6] 
for Mo isotopes, we can see the accumulated GT strength 
plotted up to the corresponding Qp energy of each iso- 
tope, which is the relevant energy range for the calcula- 
tion of the half-lives. Also shown by vertical dashed lines 
are the S n energies when they are lower than Qp. In this 
magnified scale one can appreciate the sensitivity of these 



FIG. 4: (Color online) Same as in Fig. [3l but for Mo isotopes. 



distributions to deformation and how measurements of 
the GT strength distribution from /3-decay can be a tool 
to get information about this deformation, as it was car- 
ried out in Refs. [5^, [5t|. The accumulated strength 
from the oblate shapes is in general larger than the cor- 
responding prolate profiles. The spherical distributions 
have distinct characteristics showing always as a strong 
peak at an excitation energy around 2 MeV. The pro- 
files from different shapes could be easily distinguished 
experimentally from each other. This is specially true in 
the case of the lighter isotopes ioo-i04 Zr and io4-ios MO) 
where the differences are enhanced. These isotopes are in 
principle easier to measure since they are the less exotic. 

Experimental information on GT strength distribu- 
tions in these isotopes is only available in the energy 
range below 1 MeV for the isotopes 106 - 108 Mo [goj], 110 Mo 
U, and ioo,i02,io4 Zr jg| These data can be seen in 

Figs. [5] and HI together with the QRPA calculations. 
Unfortunately, the energy region is still very narrow and 
represents only a small fraction of the GT strength rel- 
evant for the half-life determination. Clearly, more ex- 
perimental information is needed to get insight into the 
nuclear structure of these isotopes. 
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FIG. 5: (Color online) QRPA-SLy4 accumulated GT 
strengths in Zr isotopes calculated for the various equilib- 
rium shapes. In each isotope the energy range considered 
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FIG. 6: (Color online) Same as in Fig. [5] but for Mo isotopes. 



C. Half-lives and /3-delayed neutron-emission 
probabilities 

The calculation of the half-lives in Eq. (|9]) involves 
the knowledge of the GT strength distribution and of 
the Qp values. The calculation of the probability for (3- 
delayed neutron emission P n in Eq. ()13[) involves also 
the knowledge of the S„ energies. We use experimental 
values for Qp and S n , which are taken from Ref. [33| or 
from the Jyvaskyla mass database [63 |. when available. 
But in those cases where experimental masses are not 
available, one has to rely on theoretical predictions for 
them. There are a large number of mass formulas in the 



market obtained from different approaches. 

The strategy used in this work is first of all to compare 
with experiment the predictions of some representative 
mass formulas in the mass region where data are avail- 
able. According to their success in reproducing the Qp 
and S n energies, we finally adopt the most convenient 
mass formula for extrapolations to the unknown regions. 
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FIG. 7: (Color online) Experimental Qp and S n energies com- 
pared to the predictions of various mass models. 

In Fig. [7] we show this comparison for six frequently 
used mass formulas. We use the model of P. Moller et 
al. [lj| that belongs to a microscopic-macroscopic type 
of calculation. It contains a FRDM corrected by micro- 
scopic effects obtained from a deformed single-particle 
model based on a folded Yukawa potential including pair- 
ing in the Lipkin-Nogami approach. Then we use the 
ETFSI model [11] , which adopts a semi-classical ap- 
proximation to the Hartree-Fock method including full 
Strutinsky shell corrections and BCS pairing correlations. 
The label SLy4 stands for the masses calculated from 
the Skyrme force SLy4 with a zero-range pairing force 
and Lipkin-Nogami obtained from the code HFBTHO 
[6oT | and compiled in (3||. The results under the label 
Gogny have been obtained from the HFB calculations 
with the finite-range Gogny-DlS force The HFB-17 
model is one of the most recent versions of the Skyrme 
HFB mass formulas introduced by the Brussels-Montreal 
group [H,[63]. As in the case of the previous cases, SLy4 
and Gogny, this is a fully microscopic approach since it is 
based on an effective two-body nucleon-nucleon interac- 
tion. The Duflo and Zuker (DZ) mass model [57j is writ- 
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ten as an effective Hamiltonian which contains two parts, 
a monopole term and a multipole term. The monopole 
calculations are purely HF-type based on single-particle 
properties while the multipole term acts as a residual 
interaction and the calculation goes beyond HF. Its pre- 
dictive power has been recently checked [68| with a num- 
ber of tests probing its ability to extrapolate with very 
good results. In this work we use the 10-parameter ver- 
sion of the mass formula (69j . which is a simplification 
of the more sophisticated 28-parameter mass formula in 
Ref. H3j. 

In the upper panels of Fi g. [7] we can see the experimen- 
tal Qp values (black dots) (33l[63|. available for the iso- 
topes 100,102,104 Zr 104,106,108,110 Ma They ^ cQm _ 

pared with the predictions of the various mass models 
discussed above. In the lower panels we have the neutron 
separation energies S n corresponding to the daughter iso- 
topes of Nb and Tc, where we compare the measured 
energies (black dots) with the predictions of the DZ for- 
mula and SLy4 force. We have selected for consistency 
the SLy4 predictions, but also the DZ mass formula as 
one of the the most suited formula in this particular mass 
region. They agree pretty well with the measured values 
for both Qp and S n values. In what follows the results 
for half-lives and P n for the i°6,i08,iio Zr and U2,iW Mo 
will be obtained by using Qp and S n from SLy4 and DZ 
mass formula. 

In Figs. [5] and |H] we can see the dependence of the 
half-lives T X j 2 and P n values with the quadrupole de- 
formation /3. Solid lines in the lighter isotopes ( 104 Zr 
and 108 ' 110 Mo) correspond to QRPA-SLy4 calculations 
using experimental Qp and S n . In the heavier isotopes 
( 106 ~ 110 Zr and 112 ' 114 Mo), where there are no data for 
Qp and S n , solid (dashed) lines correspond to QRPA- 
SLy4 calculations using Qp and S n from SLy4 (DZ). Ex- 
perimental data are shown by horizontal dashed lines, 
where the shaded region in between corresponds to a 1- 
<j confidence level. The vertical dashed lines show the 
self-consistent quadrupole deformations for which we ob- 
tained the equilibrium shape configurations (see Figs. [T] 
and The first evidence to mention is that a spherical 
approach to these nuclei is far from the measured data, 
demanding a deformed treatment. 

In Fig. [5]we show the results for the isotopes 104_110 Zr. 
In the cases of 104 < 106 Zr we reproduce the experimental 
half-lives with oblate and prolate deformations very close 
to the self-consistent ones. In the oblate case the calcula- 
tion gives half- lives lower than experiment, while the self- 
consistent prolate deformation produces somewhat larger 
ones. Thus, the experiment would be reproduced either 
by nuclear deformations which do not produce shapes at 
equilibrium (/? « ±0.3) or by a mixing of the equilibrium 
deformations. Interestingly, similar results were obtained 
in Ref. [l4[ for 104 Zr, from the analysis of the measured 
/3-decay properties of 104 Y in terms of quadrupole de- 
formation 62 of the daughter 104 Zr (see [19| for a formal 
definition of the parameters £2 and (3). In that case, the 
quadrupole deformation needed to reproduce the mea- 
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FIG. 8: (Color online) QRPA-SLy4 /3-decay half-lives and P„ 
values for Zr isotopes as a function of the quadrupole defor- 
mation f3 compared to experiment (shaded area) . See text for 
more details. 
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FIG. 9: (Color online) Same as in Fig. [8j but for Mo isotopes. 



sured half-life and P n value was | €2 1 « 0.25, although 
oblate deformations were ruled out at that time. An im- 
portant difference between the results shown in Fig. [8] 
and those discussed in Ref. [HI is the abrupt increase 
seen in the latter for the Ti/2 and P n value for a near- 
spherical 104 Zr. These large values were mostly produced 
by the location of the GT-populated 7rgg/ 2 ( 8)J y g7/2 level at 
rather high energies in the spherical daughter 104 Zr. This 
discrepancy emphasizes the sensitivity of and P n to 
the structure details of the mother/daughter nuclei. The 
experimental P n values are only upper limits, although 
they are much larger than the typical values obtained 
theoretically. In the heavier isotopes there are no data 
and these results are then useful to see the sensitivity 
to deformation of the predictions. The spherical minima 
in the heavier isotopes predict half-lives and P n values 
much lower than the corresponding values for deformed 
shapes. In Fig. [§] we have the results for the isotopes 
i08-ii4 Mo In the cage of i08 Mo the half-life is repro- 
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duced with the self-consistent oblate deformation, while 
the prolate one generates too high half-lives. In the case 
of 110 Mo the measured half-life is well reproduced with 
both oblate and prolate equilibrium deformations. In the 
case of 108 Mo the P n value is zero since experimentally 
63] S n > Qp. For 110 Mo the P n value is not reached by 
the calculations. As in the case of the heavier Zr isotopes, 
the heavier Mo isotopes show that the half-lives for the 
spherical minima are much smaller than the correspond- 
ing half-lives for the self-consistent oblate and prolate 
shapes. In general we observe that the half-lives (P n val- 
ues) in the heavier Zr and Mo isotopes calculated with 
Qp and S n from SLy4 (solid lines) are lower (larger) than 
the results calculated with Qp and S n from DZ (dashed 
lines). 

In Figs. [TU1 and [TT1 we compare the measured /3-decay 
half-lives (upper panels) and P n values (lower panels) 
with the theoretical results obtained with the oblate, pro- 
late, and spherical equilibrium shapes. In the 10 °- 104 Zr 
isotopes we use experimental Qp and S n values, while in 
the heavier 106 - 110 Zr isotopes we use SLy4 in Fig. QJJ 
and the DZ mass formula in Fig. 111! Similarly, for the 
104_110 Mo isotopes we use experimental Qp and S n val- 
ues, while for 112 ~ 114 Mo we use SLy4 in Fig. [TUJ and 
the DZ mass formula in Fig. [TT] In the case of Zr iso- 
topes, we can see that the experimental half-life is close to 
the oblate result in 100 Zr and appear systematically be- 
tween the prolate and oblate calculations in the isotopes 
102 ' 104 ' 106 Zr. One wonders whether such result could be 
explained by the coexistence of a highly-deformed prolate 
ground-state configuration with a moderately deformed 
minimum similar to that found in 100 Zr (5ll . l52j | and (more 
speculatively) in 102 Zr [53(. The results seem to indicate 
that such weakly-deformed intruder configurations may 
have an oblate character. In the heavier isotopes 108 > 110 Zr 
the predictions of both oblate and prolate are very close 
to each other and much larger than the result obtained 
from spherical shapes. Measuring these half-lives and P n 
values will be a good opportunity to check the role of 
spherical configurations in these exotic nuclei, since the 
spherical components will lower the half-lives and P n val- 
ues by factors about 5 and 15-50, respectively. 

In the case of Mo isotopes, the experimental half-lives 
in i04,i06,i08,iio Mo tend t0 favor the blate theoretical 

results (which are indeed the ground states) over the pro- 
late ones. In the heavier 112,114 Mo isotopes, as in the case 
of the heavier Zr isotopes, oblate and prolate results are 
very similar and much larger than the spherical predic- 
tions, offering again a sensitive test to analyze the de- 
formation of these heavy nuclei, for which spectroscopic 
measurements are more difficult. Experimental P n values 
are only upper limits except for the case 110 Mo, which is 
much larger than the calculations. This implies that the 
relative GT strength contained in the energy region be- 
low S n is overestimated theoretically and therefore the 
relative contribution coming from the strength above S n 
is too small. This can be seen in Fig. [5]for 110 Mo, where 
the accumulated strength is practically flat above S n . 



The half- lives and P n values of the A ~ 110 nuclei, pre- 
dicted here for spherical configurations, would have clear 
consequences in the calculation of r-process abundances. 
In particular, the abrupt reduction of the P n values may 
contribute to fill the artificial trough around A = 110 pre- 
dicted by current r-process nucleosynthesis models. Fur- 
thermore, the confirmation of spherical shapes in these 
nuclei may be an indirect signature of the N = 82 shell 
quenching since both phenomena are predicted by the 
SLy4 force used in our calculations. 
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FIG. 10: (Color online) Measured /3-decay half-lives and P n 
values for Zr and Mo isotopes compared to theoretical QRPA 
results calculated from different shape configurations, using 
SLy4 to compute Qp and S n in the heavier isotopes. 



IV. CONCLUSIONS 

In this paper we have studied the /3-decay properties of 
neutron-deficient Zr and Mo isotopes within a deformed 
QRPA approach based on mean fields generated from 
self-consistent Skyrme Hartree-Fock calculations. In par- 
ticular, we have analyzed the experimental information 
on the half-lives and /3-delayed neutron-emission proba- 
bilities in the neutron-rich 10 °- 110 Zr and 104 ~ 114 Mo iso- 
topes in terms of the nuclear deformation. 

We have shown that the measured half-lives in Zr iso- 
topes are placed between the results obtained from the 
oblate and prolate coexistent shapes that appear very 
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FIG. 11: (Color online) Same as in Fig. \W\ but using DZ 
instead of SLy4 to compute Qp and S n in the heavier isotopes. 



in the heavier isotopes 112 ' 114 Mo, the predicted half-lives 
for both shapes are very close and larger than the spher- 
ical ones. Thus, comparison with experimental half-lives 
indicates that in some cases (Mo isotopes) a single shape 
accounts for this information, while in other cases (most 
of the Zr isotopes) a more demanding treatment in terms 
of mixing of different shapes seems to be more appropri- 
ate. P n values are in general not well reproduced, al- 
though experimentally only upper limits are measured 
in most cases. Hence, it will be certainly worth mea- 
suring those heavier isotopes and check whether they are 
properly described by the deformed shapes, or a spherical 
component is needed as well. 

Nevertheless, one should keep in mind that half-lives 
are integral properties that collect all the information of 
the decay in a single number and does not tell us about 
the detailed internal structure of the GT strength dis- 
tribution, much more sensitive to the nuclear structure. 
From a more detailed analysis of the GT strength dis- 
tributions in the Qp energy range accessible in /3-decay, 
we have shown that the differences between the predic- 
tions of the different nuclear shapes could be clearly dis- 
tinguished experimentally. Although these spectroscopic 
measurements are at present not feasible because of the 
still low production rates of exotic nuclei at modern ra- 
dioactive beam facilities, they will provide in the future 
precise tests of the nuclear structure in exotic nuclei. 



close in energy in the PECs. The predicted half-lives in 
the heavier Zr isotopes 108 > 110 Zr, where there are no ex- 
perimental data yet, are however very close to each other 
for both oblate and prolate shapes and much larger than 
the predictions from the spherical shapes. On the other 
hand, the measured half-lives in Mo isotopes agree bet- 
ter with the calculations from oblate shapes, which are 
lower than the corresponding prolate ones. Once more, 
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